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Magnetoresistive effects are usually invariant upon inversion of the magne¬ 
tization direction. In noncentrosymmetric conductors, however, nonlinear 
resistive terms can give rise to a current dependence that is quadratic in 
the applied voltage and linear in the magnetization. Here we demonstrate 
that such conditions are realized in simple bilayer metal films where the 
spin-orbit interaction and spin-dependent scattering couple the current- 
induced spin accumulation to the electrical conductivity. We show that 
the longitudinal resistance of Ta|Co and Pt|Co bilayers changes when re¬ 
versing the polarity of the current or the sign of the magnetization. This 
unidirectional magnetoresistance scales linearly with current density and 
has opposite sign in Ta and Pt, which we associate to the modification of 
the interface scattering potential induced by the spin Hall effect in these 
materials. Our results suggest a new route to detect magnetization switch¬ 
ing in spintronic devices using a two-terminal geometry, which applies also 
to heterostructures including topological insulators. 
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The effects of the magnetization on the electric conductivity of metals have been studied 
for a long time 1 , providing understanding of fundamental phenomena associated to elec¬ 
tron transport and magnetism as well as manyfold applications in sensor technology. The 
anisotropic magnetoresistance (AMR), the change of the resistance of a material upon ro¬ 
tation of the magnetization, is a prominent manifestation of spin-orbit coupling and spin- 
dependent conductivity in bulk ferromagnets 2,3 . In thin film heterostructures, the additional 
possibility of orienting the magnetization of stacked ferromagnetic layers parallel or antipar¬ 
allel to each other gives rise to the celebrated giant magnetoresistance (GMR) effect 4,5 , which 
has played a major role in all modern developments of spintronics. Together with the early 
spin injection experiments 6,7 , the study of GMR revealed how nonequilibrium spin accu¬ 
mulation at the interface between ferromagnetic (FM) and normal metal (NM) conductors 
governs the propagation of spin currents 8 ” 11 and, ultimately, the conductivity of multilayer 
systems 10,12 . 

Recently, it has been shown that significant spin accumulation at a FM|NM interface can 
be achieved using a current-in-plane (CIP) geometry owing to the spin Hall effect (SHE) 
in the NM 13 . When FM is a metal and NM is a heavy element such as Pt or Ta, the spin 
accumulation is strong enough to induce magnetization reversal of nm-thick FM layers at 
current densities of the order of j — 10' —10 8 A/cm 2 (Refs. 14-16). When FM is an insulator, 
such as yttrium iron garnet, the SHE causes an unusual magnetoresistance associated to the 
back-flow of a spin current into the NM when the spin accumulation fj, s ~ (j x z) is aligned 
with the magnetization of the FM, which increases the conductivity of the NM due to 
the inverse SHE 1 '” 19,25 . This so-called spin Hall magnetoresistance (SMR) is characterized 
by R v < R z R x , where R 1 is the resistance measured when the magnetization (M) is 
saturated parallel to i — x,y,z, and differs from the conventional AMR in poly crystalline 
samples, for which R — R x — ( R y,z — R X )[M • j ] 2 and R v « R z < R x (Ref. 3). 

In this work, we report on a new magnetoresistance effect occurring in FM|NM bilayers 
with the NM possessing large SHE. The effect combines features that are typical of the 
current-in-plane (CIP) GMR and SHE, whereby the spin accumulation induced by the SHE 
in the NM replaces one of the FM polarizers of a typical GMR device. Differently from GMR, 
however, this effect introduces a nonlinear dependence of the resistance on the current, which 
gives it unique unidirectional properties: the resistivity changes when reversing either the 
sign of the magnetization or the polarity of the current, increasing (decreasing) when the 
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SHE-induced non equilibrium magnetization at the FM|NM interface is oriented parallel 
(antiparallel) to the magnetization of the FM, as illustrated in Fig. la,b. We associate 
this phenomenon to the modulation of the FM|NM interface resistance due to the SHE- 
induced spin accumulation, which gives rise to a nonlinear contribution to the longitudinal 
conductivity that scales proportionally with the current density and has opposite sign in 
Pt and Ta. Contrary to the linear magnetoresistive effects, including the AMR, GMR, and 
SMR described above, which are even with respect to the inversion of either the current or 
magnetization owing to the time reversal symmetry embodied in the Onsager’s reciprocity 
relations, here we observe R( j. M) = — R(— j, M) = —R( j, — M), providing a unidirectional 
contribution to the magnetoresistance in simple bilayer systems. 

Sample layout 

The samples studied in this work are Pt(l-9nm)|Co(2.5nm) and Ta(l-9nm)|Co(2.5nm) films 
with spontaneous in-plane magnetization, capped by 2 nm AlO^ and patterned in the shape 
of Hall bars of nominal length l — 20 — 50 gm, width w — 4 — 10 /ini, and l/w — 4, as 
shown in Fig. lc. Additional control experiments were carried out on single Co, Ta, and Pt 
films, and Ta(l,6nm)|Cu(2,4,6nm)|Co(2.5nm) trilayers. To characterize the magnetic and 
electrical properties of these layers we performed harmonic measurements of the longitudinal 
resistance (R, see Supplementary Information) and transverse Hall resistance (Rh) 16,2022 as 
a function of a vector magnetic field defined by the polar and azimuthal coordinates 6b and 
ipB- The measurements were carried out at room temperature by injecting an ac current of 
frequency co/2n — 10 Hz and simultaneously recording the first (iC) and second harmonic 
resistance (R^) between the contacts shown in Fig. lc while rotating the sample in a uniform 
magnetic field of 1.7 T. Here, R w represents the linear response of the samples to the current, 
that is, the conventional resistance. In order to include the different magnetoresistive angular 
dependencies in a single expression we write this term as 

Ru — R z + ( R x — R z ) sin 2 6 cos 2 </> + ( R y — R z ) sin 2 6 sin 2 ip , (1) 

where 6 and are the polar and azimuthal angles of M, as schematized in Fig. 2a. i? 2 w de¬ 
scribes resistance contributions that vary quadratically with the applied voltage and includes 
the current-induced changes of resistivity that are the main focus of this work. 

Magnetoresistance measurements 
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Figure 2b and c show the resistance of Ta(6nm)|Co(2.5nm) and Pt(6nm)|Co(2.5nm) lay¬ 
ers during rotation of the applied field in the xy, zx, and zy planes. We observe a size¬ 
able magnetoresistance (MR) in all three orthogonal planes and R x > R z > R y for both 
samples, in agreement with previous measurements on Pt|Co films 23 ’ 24 . The resistivity 
of Ta(6nm)|Co(2.5nm) [Pt(6nm)|Co(2.5nm)] is 108.9 (36.8) //Mem and the MR ratios are 
(R x — R z )/R z = 0.09% [0.05%] and (R z — R v )/R z — 0.12% [0.53%], showing a large SMR-like 
behavior compared to Ta|YIG and Pt|YIG 17,25 . The solid lines represent fits to the MR using 
Eq. 1 and 9 simultaneously measured via the anomalous Hall resistance (see Supplementary 
Information). 

In addition to the linear resistance, we measure an unexpected nonlinear resistance, R 2w , 
which has a different angular dependence compared to R w and opposite sign in Pt and Ta, 
as shown in Fig. 2d and e. By fitting the curves with respect to the angles 9 and ip (solid 
lines), we find that R 2uj ~ siriO sin p ~ M y . In the following, we discuss the type of nonlinear 
effects that can give rise to such a symmetry. 

Spin-orbit torques and thermoelectric contributions 

First, we consider oscillations of the magnetization due to the current-induced spin-orbit 
torques (SOT) 16,20-22,26 . As the SOT are proportional to the current, ac oscillations of the 
magnetization can introduce second-order contributions to R 2u due to the first order MR 
described by Eq. 1. However, as shown in detail in the Supplementary Information, the 
SOT-induced signal is not compatible with the angular dependence of R 2uj . Both the field¬ 
like and antidamping-like SOT (as well as the torque due to the Oersted held) vanish for 
M || y, where |i? 2 w| is maximum. Moreover, the field-like SOT is small in 2.5 nm thick Co 
layers 22 , whereas the antidamping SOT can only induce variations of R 2co in the zx plane 
with maxima and minima close to 9b = 0° and 180°, which we observe to be small and more 
pronounced in Pt|Co relative to Ta|Co (Fig. 2d and e bottom panels). 

Second, we analyze the influence of thermal gradients (VT) and related thermoelectric 
effects. The anomalous Nernst (ANE) and spin Seebeck effect (SSE) 2 '’ 28 , both inducing 
a longitudinal voltage proportional to y 2 (M x VT), can give rise to a similar angular de¬ 
pendence as observed for R 2uj when VT || z (see Supplementary Information). Here, we 
find that thermoelectric voltages are negligible in Pt|Co, in agreement with the very small 
thermal gradients reported for this system 22 . In Ta|Co, on the other hand, the much larger 
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resistivity of Ta relative to Co results in a higher current flowing in the Co layer and a 
positive VT. In such a case, the second harmonic signal of thermal origin, R%J, can be 
simply estimated from its transverse (Hall) counterpart scaled by the geometric factor l/w 
when the magnetization is tilted in the x direction, and subtracted from the raw signal. 
Accordingly, we find that Rjj — 5 mO in Ta(6nm)|Co(2.5nm), which accounts for only 
about 50% of the total R 2bJ reported in Fig. 2. The same procedure applied to Pt|Co gives 
Rjj mil of the order of 5% of the total R 2uj , whereas in the control samples lacking a heavy 
metal we find uniquely a signal of thermal (ANE) origin. We conclude that there is an 
additional magnetoresistive effect in the Pt|Co and Ta|Co bilayers that cannot be accounted 
for by either current-induced magnetization dynamics or thermoelectric voltages. 

Unidirectional spin Hall magnetoresistance 

The symmetry as well as the opposite sign of the nonlinear resistance in Ta|Co and Pt|Co 
suggest that it is related to the scalar product of the magnetization with the SHE-induced 
spin accumulation at the FM|NM interface, (j xz) -M, giving rise to a chiral MR contribution 
rusmr ^ j x M. This relation describes the general features expected from a unidirectional 
magnetoresistance driven by the spin Hall effect (USMR). We note that this MR contribution 
depends on the current direction and that the resistance of the bilayer increases when the 
direction of the majority spins in the FM and the spin accumulation vector are parallel to 
each other, and decreases when they are antiparallel. This may appear counterintuitive at 
first sight, considering that the conductivity of Co is larger for the majority spins. However, 
as we will discuss later, this behavior is consistent with the theory of GMR in FM|NM|FM 
heterostructures 10 ’ 29 ’ 30 when only a single FM|NM interface is retained and the SHE is taken 
into account. 

To investigate further the USMR we have measured R 2u as a function of an external 
magnetic held applied parallel to y and current amplitude. Figure 3a shows that R 2u > is 
constant as a a function of held for Ta|Co as well as for Pt|Co (Fig. 3b) and reverses 
sign upons switching of the magnetization from the y to the — y direction. In the Pt|Co 
case we observe also two spikes, which we attribute to the magnetization breaking into 
domains at low held and giving rise to dynamic effects on the domain walls 16 . Note that the 
magnetization of Pt|Co is not fully saturated below 0.65 T due to the large perpendicular 
magnetic anisotropy of this system, differently from Ta|Co (Supplementary Information). 
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Figure 3c shows the current dependence of R. 2 RMR — R 2 u ~ R%J (R- 2 u MR ~ R- 2 u for Pt|Co) 
obtained by taking the average of the data measured at fields larger than | ± 1| T. R\ ,RMR 
is linear with the injected current density and converges to zero within the error bar of the 
linear fit (black lines). 

To verify the role of the interfacial spin accumulation due to the SHE we examined the 
dependence of the USMR on the thickness of the NM. Figure 4a and b show the absolute 
change of sheet resistance A R usmr = R^ rmr (± M, ±j) — R rrmr (±M. =pj) and the relative 
change of resistivity /\r usmr /R measured at constant current density as a function of the 
Ta and Pt thickness. Both curves exhibit qualitatively similar behavior: an initial sharp 
increase below 2-3 nm and a gradual decrease as the NM layer becomes thicker. We note 
that the USMR signal is almost absent in Ta(lnm)|Co, contrary to Pt(lnm)|Co, which we 
attribute to the oxidation of Ta when deposited on SiCU and its consequent poor contribution 
to electrical conduction. The initial increase of the USMR is consistent with the increment 
of the spin accumulation at the FM|NM interface as the thickness of the NM becomes 
larger than the spin diffusion length, which is of the order of 1.5 nm in both Ta and Pt 25,31 . 
Moreover, we observe that the decline of the signal in the thicker samples is stronger in Pt|Co 
than in Ta| Co. This behavior is consistent with Pt gradually dominating the conduction 
due to its low resistivity, and a smaller proportion of the current experiencing interface 
scattering in Pt|Co. Conversely, the high resistivity of Ta shunts the current towards the 
Co side, increasing the relative proportion of the current affected by scattering at the Ta|Co 
interface. 

As an additional check to validate these arguments we have performed measurements on 
single Ta(6nm), Pt(6nm), and Co(8nm) layers as well as on Ta(l,6nm)|Cu(2,4,6nm)|Co(2.5nm) 
trilayers, all capped by 2 nm A10 x . The USMR is absent in the Ta, Pt, and Co single layers, 
which also excludes any self induced magnetoresistive effect 32 and proves the essential role 
of the FM|NM interface. On the other hand, we find a sizable USMR when a Cu spacer is 
inserted between Ta and Co, which excludes proximity-induced magnetization as a possible 
cause for the USMR (see Supplementary Information). 

Discussion 

Based on the analysis presented above, we conclude that the current-induced spin accumu¬ 
lation creates an additional spin-dependent interface resistance that adds or subtracts to the 
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total resistance depending on sign of the cross product j x M. Given the in-plane geometry, 
the interpretation of this effect requires a Boltzmann equation approach to model the spin- 
and momentum-dependent reflection and transmission of electrons at the FM|NM interface, 
equivalent to extending the theory of CIP-GMR 29,30 beyond first order and including the 
SHE. However, a qualitative understanding of the USMR can be derived also by introducing 
a two current series resistor model and an interface resistance term proportional to the SHE- 
induced shift of the electrochemical potential between the FM and NM layers. The latter 
can be calculated using a one-dimensional drift-diffusion approach 9,10,17 . We consider two 
separate conduction channels for the minority (spin f) and majority (spin J,) electrons. As in 
bulk FM, scattering at the interface is spin-dependent due to the unequal density of majority 
and minority states near the Fermi level, which, in most cases, leads to a larger resistance 
for minority electrons relative to majority electrons: r 7 > r 7 . This resistance difference is 
at the heart of the GMR effect, both in the CIP 29,30,33 and the current-perpendicular-to- 
plane (CPP) geometry 10,34 . Additionally, when an electric current flows from a FM to a NM 
or viceversa, another resistance term appears due to the conductivity mismatch between 
majority and minority electrons on opposite sides of the junction, which results in spin ac¬ 
cumulation (Refs. 8 and 9). This so-called ”spin-coupled interface resistance”, plays a role in 
CPP-GMR as well as in local and nonlocal spin injection devices'’ 10,35 , whereas in the CIP 
geometry it is normally neglected because there is no net charge flow across the interface 
and the spin accumulation is assumed to be zero. If we take the SHE into account, however, 
the transverse spin current flowing between the NM and the FM induces a splitting of the 
spin-dependent electrochemical potentials /A and p^ and a net interfacial spin accumulation 
p s = p^ — /U 7 , which is given by 


A 4 s N At 


,o 


tanh 




1 + 


n /i 


( 2 ) 


where /j° v = 2 e Pn^n Osh j is the bare spin accumulation due to the SHE that would occur 
in a single, infinitely thick NM layer, 0 SH the spin Hall angle of the NM, Pat,f and Xn,f are 
the resistivity and spin diffusion length of the NM and FM, respectively, and — (r 7 + r 7 )/4 
is the interface resistance 10 . Moreover, the same effect induces a shift Ap — Pn — Pf of the 
electrochemical potential pn,f = (p% f + A%f)/2 of the NM relative to the FM: 
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where 7 = (H" — r^)/(r^ + H") and r = W^(l — P 2 ) tanh A). Figure 5a and b show a 
graphical representation of p s and Adetails about the derivation of Eqs. 2 and 3 are 
given in the Supplementary Information. A key point is that A/pv depends on the product 
Pd SH j , as the USMR, and is linked with the spin-dependent scattering potential that gives 
rise to different transmission coefficients for majority and minority electrons at the FM|NM 
interface 36 . We can thus draw the following qualitative interpretation of the USMR: when 
the nonequilibrium magnetization induced by the SHE and the magnetization of the FM 
are parallel to each other (<5m || M), the transmission of t (i) electrons across the interface 
is reduced (enhanced) by the accumulation of majority electrons at the FM|NM boundary, 
due to the conductivity mismatch of t and \. spins in the two materials. Likewise, when 
Am || —M, the transmission of J, (t) electrons across the interface is reduced (enhanced) since 
minority electrons accumulate at the FM|NM boundary. The overall effect is a modulation of 
the interface resistance of the t and J, spin channels by a nonlinear term ±r s , as schematized 
in Fig. 5c. This two current (t and series resistor model shows that such a modulation 
leads to a resistance difference between the two configurations AR^ mr 2 r s j. 

Accordingly, using realistic values of r^, pjy, and p F for Ta|Co and Pt|Co, we fit the 
dependence of the USMR on current and NM thickness to the following phenomenologi¬ 
cal expression (see Supplementary Information): AR usmr = A tanh W^/(l + Rfi/Rn) 2 , 
where A is a parameter proportional to Pp° SN representing the amplitude of the effect, Rpi 
is the effective resistance of the FM and interface regions, and Rn = ppfl/{wtpi) is the 
resistance of the NM. The denominator accounts for the decreased fraction of electrons 
that scatter at the interface as the thickness of the NM increases. Similarly, we obtain 
A R usmr /R = ( A/R ) tanh A^/(l + R F i/Rn). As shown in Fig. 4, these simple expressions 
fit A R usmr and AR usmr /R remarkably well, providing also values of A p t = 1.1 nrn and 
\ Ta = 1.4 nm that are in agreement with previous measurements 25,31 . Our model thus 
captures the essential features of the USMR, namely its sign, angular dependence, and pro¬ 
portionality to the current. Detailed calculations including realistic scattering parameters 
within a nonlinearized Boltzmann approach including the SHE 3 ' should be able to account 
for quantitative aspects of the USMR in different materials. We stress also that the USMR 
is not uniquely linked to the SHE but may arise also due to other sources of nonequilib¬ 
rium spin accumulation, such as the Rashba effect at FM|NM interfaces and topological 
insulators 37-40 as well as the anomalous Hall effect in FM. 



Conclusions 

The existence of a nonlinear magnetoresistive term proportional to j x M has both funda¬ 
mental and practical implications. Identifying which symmetries survive the breakdown of 
the Onsager relationships in the nonlinear regime is central to the understanding of elec¬ 
tron transport phenomena, particularly in mesoscopic and magnetic conductors where such 
effects can also have thermoelectric and magneto-optical counterparts 41,42 . In this respect, 
the USMR shows that the longitudinal conductivity has an antisymmetric Hall-like com¬ 
ponent that has so far remained unnoticed. We expect such a component to be a general 
feature of noncentrosymmetric magnetic systems with strong spin-orbit coupling. We note 
also that the USMR differs from the nonlinear MR observed in chiral conductors, such as 
metal helices 42 and molecular crystals 43 , which is proportional to j • M. 

In the field of spintronics, nonlinear interactions between spin and charge are emerging as 
a tool to detect spin currents 44 and thermoelectric 45 effects, as well as magnetization reversal 
in dual spin valves 46 . Although the USMR is only a small fraction of the total resistance, its 
relative amplitude is of the same order of magnitude of the spin transresistance measured 
in nonlocal metal spin valves 7,35 , which is a method of choice for the investigation of spin 
currents. The thermoelectric counterpart of the USMR, related to the spin Nernst effect, 
may be used to detect lieat-induced spin accumulation by modulation of the magnetization 
rather than an electric current. We note that the electric field created by the USMR is of 
the order of 2 V/m per 10 7 A/cm 2 , which is comparable to the ANE 22 and three orders of 
magnitude larger than the typical electric fields due to the SSE 2 ' ,2S . 

In terms of applications, the USMR may be used to add 360° directional sensitivity to 
AMR sensors, which are widely employed for position, magnetic field, and current sensing, 
and already include built-in modulation circuitry for accurate resistance measurements. 
Most interestingly, the USMR shows that it is possible to realize two-terminal spintronic 
devices where switching is performed by SOT 14,15 and reading by a resistance measurement. 
Such a scheme involves only one FM layer and minimum patterning effort. Finally, we 
believe that there is substantial room for improving the amplitude of the USMR to levels 
closer to the AMR, either by material or heterostructure engineering. In particular, the 
USMR could increase significantly in magnetic topological insulator structures due to the 
very large spin accumulation and reduced bulk conductance reported for these systems 47,48 . 
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Methods 

Sample preparation. The Pt(l-9nm)|Co(2.5nm)|A10 x (2nm) and Ta(l-9nm)|Co(2.5nm)|A10 x (2nm) 
layers were grown by dc magnetron sputtering on thermally oxidized Si wafers. The depo¬ 
sition rates were 0.185 nm/s for Pt, 0.067 nm/s for Ta, 0.052 nm/s for Co, and 0.077 nm/s 
for Al. The deposition pressure was 2 mTorr and the power was 50 W for all targets. The 
A1 capping layers were oxidized in-situ by a 7 mTorr O 2 plasma at 10 W during 35 s. The 
layers were subsequently patterned into 6-terminals Hall bars by using standard optical 
lithography and Ar milling procedures. The Hall bar dimensions are w for the current line 
width, w/2 for the Hall branch width, and l — Aw is the distance between two Hall branches, 
where w varies between 4 and 10 p m. 

Characterization. All layers possess spontaneous isotropic in-plane magnetization. To 
determine the saturation magnetization of Co we have performed anomalous Hall effect 
measurements on an 8 nm-thick Co reference sample with B ext || z. The field required to 
fully saturate M out-of-plane is about 1.5 T, which, assuming that perpendicular magnetic 
anisotropy is negligible in this layer, is close to PqM s expected of Co. Similar measurements 
on Ta(6nm)|Co(2.5nm) and Pt(6nm)|Co(2.5nm) layers give saturation fields of 1.45 T and 
0.8 T, respectively. This is attributed to the small (large) perpendicular interface anisotropy 
contribution of the Ta|Co (Pt|Co) interface, reducing the field required to saturate the mag¬ 
netization out-of-plane. 4-point resistivity measurements on single Co(8nm), Ta(6nm) and 
Pt(6nm) layers yield pc 0 — 25.3 pQ, cm, p Ta — 237 pQcm and p Pt — 34.1 pflcm, in line with 
expectations for Pt and Co thin films, and the /3-phase of Ta. The magnetoresistance and 
Hall voltage measurements were performed at room temperature by using an ac current 
/ = /o sin cut where cu/27t = 10 Hz, generated by a Keithley-6221 constant current source. 

For the data reported in Fig. 2 the peak amplitude of the injected ac current was set to 
8.5 mA, corresponding to a nominal current density of j — 10' A/cm 2 . In other measure¬ 
ments with different device size and thickness, the current was adapted to have the same 
current density. The longitudinal and transverse voltages were recorded simultaneously by 
using a 24-bit resolution National Instruments PXI-4462 dynamic signal analyser, during 
10 s at each angle position in a uniform external field of 1.7 T. The rotation of the sample 
was provided by a motorized stage having a precision of 0.02 degrees. The acquired voltages 
were fast Fourier transformed (FFT) to extract the first and second harmonic voltage signals 
14 and V 2 u - The corresponding resistances are given by — V u /Iq and = V^/Iq (peak 
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values). 
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FIG. 1. Illustration of the unidirectional spin Hall magnetoresistance elfect and sample 
layout, a, Parallel alignment of the SHE-induced nonequilibrium magnetization at the FM|NM 
interface with the magnetization of the FM increases the resistivity of the bilayer, b, Antiparallel 
alignment decreases the resistivity. The arrows indicate the direction of the spin magnetic mo¬ 
ment. c, Scanning electron micrograph of a sample and schematics of the longitudinal resistance 
measurements. 
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FIG. 2. Linear and nonlinear magnetoresistance, a, Geometry of the measurements, b, First 
harmonic resistance of Ta(6nm)|Co(2.5nm) and c, Pt(6nm)|Co(2.5nm) measured with a current 
density of j — 10' A/cm 2 . d,e Second harmonic resistance measured simultaneously with b,c. 
The dimensions of the Hall bars are l — 50 /tin and w — 10 /urn. 
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FIG. 3. Field and current dependence of the nonlinear magnetoresistance. of 

a, Ta(6nm)|Co(2.5nm) and b, Pt(6nm)|Co(2.5nm) recorded during a field sweep along y with a 
current density of j = 1.2 • 10 7 A/cm 2 , c, Current dependence of R 2 RMR . The solid lines are 
fits to the data. The slope gives the amplitude of the USMR, which is 5.5 mQ and 1.25 rriQ per 
10 7 A/cm 2 for these samples. The thermal contribution R^J has been subtracted from the Ta|Co 
data. The dimensions of the Hall bars are l — 50 /nn and w — 10 /im. 
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FIG. 4. USMR as a function of NM thickness, a, Sheet resistance ;1/ ^ as a function of 

Ta (squares) and Pt (circles) thickness measured at constant current density j — 10' A/cm 2 . The 
Co layer is 2.5 nm thick in all samples, b, Normalized resistance A R 2 ^ mr /R. The solid lines are 
fits to the data according to the model described in the text. 
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FIG. 5. Modulation of the spin accumulation, spin-dependent electrochemical poten¬ 
tial, and interface resistance by the SHE. Profile of the electrochemical potential of majority 
(fjJ, blue lines) and minority (/rq red lines) electrons in proximity of the FM|NM interface for a, 
positive and b, negative current. The electrochemical potential of the NM shifts relative to that 
of the FM as indicated by the green arrow. The direction of the magnetization is M || y in both 
panels. In our notation, the majority spins are oriented antiparallel to M and, likewise, the non- 
equilibrium magnetization induced by the SHE has opposite sign relative to fi s . Reversing M is 
equivalent to exchanging /v and /jt and inverting A/u The sign of the SHE corresponds to Pt|Co; 
the parameters used to calculate /A and fjr are given in the Supplementary Information, c, Two 
current series resistor model of the USMR corresponding to panel a (top, higher resistance) and b 
(bottom, lower resistance). Note that the resistances r^ and H" may be generalized to include also 
the bulk spin-dependent resistances of the FM layer. 
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SI 1. HARMONIC ANALYSIS OF THE LONGITUDINAL AND HALL RESISTANCES 

The longitudinal resistance ( R ) and transverse (Hall) resistance (R H ) are measured by applying an ac current 
I = Iq sin(wt) of constant amplitude Iq and frequency w/27t = 10 Hz and recording the ac longitudinal (V) and 
transverse ( V H ) voltage, respectively. Ohm’s law for a current-dependent resistance reads V(I) = R(I) ■ /osin(aA). 
Assuming that the current-induced resistance changes are small with respect to the linear resistance of the sample, 
R( 0), we can expand R and V as follows: 

R(I) = R( 0) + ^dl = R( 0) + i/ofsinM , (SI) 

and 

1 dB 

V (/) = IoR{0) sin(u;f) + - Iq — sin(2wf), (S2) 

where the longitudinal voltage consists of first and second harmonic terms that scale with Iq and Iq , respectively. 
Accordingly, we define the first and second harmonic resistance as R u = R( 0) and R 2lj = Analogue expansions 

apply to the first and second harmonic Hall resistances, R% and R 2bJ . The first harmonic terms can be written as 

Ru = R z + ( R x — R z ) sin 2 8 cos 2 ip + (R y — R z ) sin 2 9 sin 2 ip , (S3) 

R„ = Rahe cos 8 + Rphe sin 2 6 sin(2<p), (S4) 

where 9 and ip are the polar and azimuthal magnetization angles, R l is the longitudinal resistance measured when 
the magnetization is saturated parallel to the direction i = x, y, z , and Rahe and Rphe are the anomalous Hall and 
planar Hall coefficients, respectively. Note that Eqs. S3 and S4 simply represent the conventional magnetoresistance 
(MR) and Hall resistance. The second harmonic terms, on the other hand, include all the contributions to V and 
V H that vary quadratically with the current, specifically: the MR and Hall resistance changes due to current-induced 
spin-orbit torques (SOT) and Oersted field si S4 , thermoelectric effects S5,S6 , and nonlinear resistive terms. We have 
previously studied the influence of the SOT and thermal gradients (VT) on the transverse voltage, showing that 
the corresponding second harmonic Hall resistances, R^ SOT and R^) VT , can be separately measured due to their 
different symmetry and field dependence 36 . In this paper we show there is an additional nonlinear contribution that 
adds to the longitudinal resistance, which we call unidirectional spin Hall magnetoresistance (USMR). We thus have 


B „ r / JjSOT , nVT , nUSMR 

H2ui oc 1o(K 2uj + U 2ui + n 2u] 


), 


(S5) 


R 


H 

2oj 


oc Iq(R 


H,SOT 

2uj 


+ R 


H,VT 

2uj 


(S6) 


The SOT terms are due to the current-induced oscillations of the magnetization that modulate the MR and Hall 
resistance through the dependence of the angles 9 and ip on the current. These terms read, for the longitudinal and 
transverse case, as 


, d sin 2 6 


d sin 2 ip 


= Io[(R? - R z ) cos 2 ip + {R v - R z ) sin 2 tp - + I 0 [{R y - R z ) sin 2 8 - (R x - R z ) sin 2 8]^—^, (S7) 


dl 


dl 


II 


H,SOT 
2 uj 


= Iq(Rahe — 2Rphe cos 9 sin 2<p) 


dcos9 

dl 


+ IqRphe sin 2 9 


d sin 2 ip 
dl 


(S8) 


The thermoelectric terms are due to Joule heating and corresponding quadratic increase of the sample temperature 
with current, which gives rise to temperature gradients 36 : 


VT oc Iq sin 2 (wt)i?o = ~ cos(2wt)]i?o • 


(S9) 
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Such thermal gradients can give rise to the anomalous Nernst S7,S8 , spin Seebeck S9 and magneto-thermopower effects S7 
with distinct angular dependencies. Their functional form in the longitudinal and transverse geometries is given in 
Sect. SI 4 for an arbitrary VT. 

Finally, the modulation of the ferromagnet/normal metal (FM/NM) interface resistance by the spin Hall effect 
(SHE) gives rise to a second harmonic term 

R^ mr cx I 0 e S HM y , (S10) 

where 9sh is the spin Hall angle of the NM and M y the y component of the magnetization of the FM (see Figure 1 
of the main text for the definition of the spatial coordinates). 


SI 2. DETERMINATION OF THE POLAR MAGNETIZATION ANGLE USING THE HALL 

RESISTANCE 


If the magnetization is fully aligned to the external field, the polar and azimuthal angles of the M and B vectors 
coincide (9 = 9 b and ip = <pb)- Even for magnetic fields of the order of a few Tesla, however, the demagnetizing 
field induces a significant deviation of 9 from 9b , while, in the case of easy plane anisotropy of interest here, <p = ips 
remains true. In order to fit the angular dependence of the MR in the zx, zy, and xy planes using Eq. S3, we therefore 
need an independent measurement of 8. This is easily performed by measuring the Hall resistance as a function of 8b 
in the zx (ip = 0) or zy (<p = 7r/2) planes, as R% simplifies to R% = Rahe cos 6 in this case. The saturation value of 
the anomalous Hall resistance, Rahe , is obtained as shown in Fig. SI, where we report the Hall resistance for the zy 
angular scans corresponding to (a) Ta(6nm)|Co(2.5nm) and (b) Pt(6nm)|Co(2.5nm), simultaneously recorded during 
the magnetoresistance (MR) measurements shown in Fig. 2 of the main text. A nearly identical R^ is detected in the 
zx plane for both samples. The insets show measured as a function of external field parallel to z, showing that 
the magnetization is technically saturated at a field of ±1.7 T. Note that the saturation field of Ta|Co (~ 1.45 T) is 
larger than that of Pt|Co (~ 0.8 T) because of the strong perpendicular magnetic anisotropy (PMA) of Pt|Co, which 
counters the demagnetizing field of the Co layer. 


a Ta(6)/Co(2.5)/AIOx b Pt(6)/Co(2.5)/AIOx 




FIG. SI. Anomalous Hall resistance measurements. a, = RahecosO of Ta(6nm)|Co(2.5nm) and b, 

Pt(6nm)|Co(2.5nm) measured during a zy scan at fixed external field B = 1.7 T. A small constant offset due to misalign¬ 
ment of the Hall branches has been subtracted from both curves. Insets: R x f as a function of out-of-plane field. The features 
around B ex ± = 0 T are due to planar Hall effect contributions. 


SI 3. INFLUENCE OF SPIN-ORBIT TORQUES ON THE SECOND HARMONIC MR 

MEASUREMENTS 

An electric current flowing in the plane of a FM/NM bilayer generates two qualitatively different types of SOT: 
a field-like (FL) torque tfl ~ M x y, and an antidamping-like (AD) torque tad ~ M x (y x M), where y is the 
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FIG. S2. Macrospin simulations of and RY rmr - a, xy scan, b, zx scan, and c, zy scan simulations. The simulated 

planar Nernst signal is also shown for comparison in b. Note that the SOT signals are approximately a factor 20 smaller 
compared to the USMR in the xy and zy plane. This is due to the fact that the magnetization becomes less susceptible to the 
current-induced torques when tb is much larger compared to tad and tfl- 


in-plane axis perpendicular to the current flow direction x (Refs. SI, S4, S10-S12). The action of these torques is 
equivalent to that of two effective fields Bfi~Mx tfl and B ad ~ M x tad perpendicular to the instantaneous 
direction of the magnetization . 

Most of the studies using harmonic measurements so far have employed the Hall effect to reveal and quantify the 
S0T S1 ~ S4,S13 . However, SOT also give rise to second harmonic MR signals, as noted in Ref. S3. Here we 

analyze the angular dependence of based on symmetry arguments and macrospin simulations of the MR. To 

perform the simulations we compute the equilibrium position of M by considering the sum of all the torques acting 
on it: the torque due to the external field (tb), PMA (tpma), demagnetizing field (tdem), as well as the field-like 
SOT, including the torque due to the Oersted field, ( tfl ), and the antidamping-like SOT (tad)- At equilibrium, 
the position of the magnetization is defined by: 


Tb +T PMA + T DEM + T F L + T AD = 0 (Sll) 

We numerically solve the above equation for the equilibrium angles of the magnetization 9 and ^asa function of 
external field and current, using measured values of the following parameters for Ta(6nm)|Co(2.5nm): R x — R v = 
1.36 Q, R z - R y = 0.82 ft, R H = 1 Q, R PNE = 0.2 mf2, t F l = 0.24 mT, t A d = 0.6 mT, t DEM = 1-45 T, 
and t pma ~ 0. Figures S2a-c compare the simulations of the second harmonic MR due to the SOT (Eq. S7) and 
USMR (Eq. S10) in the xy 1 zx, and zy planes for an external field B ext = 1.7 T and current density j = 10 7 A/cm 2 . 
The simulations show that the FL torque can give rise to a second harmonic MR contribution in the xy and zy 
scans whereas the contribution due to the AD torque is equally zero in both cases. In the zx scan we have the 
opposite behavior, i.e., that the FL torque gives no contribution, whereas the AD torque gives rise to a non-zero 
signal. However, in all cases the MR contributions with SOT origin have a distinct symmetry with respect to the 
measurements reported in Fig. 2 of the main text. Moreover, = 0 when either ips or 9b = 7r/2 excluding the 

influence of SOT on the USMR measurements. 


SI 4. THERMOELECTRIC EFFECTS 

Thermoelectric effects must be carefully considered in the measurement of electrical signals due to unavoidable 
Joule heating and consequent temperature gradients in the sample. The thermal gradient can have an arbitrary 
direction depending on the device geometry, stacking order of the layers, inhomogeneous current flow due to resistivity 
inhomogeneities in each device and give rise to a second harmonic signal, RYJ", due to their quadratic dependence on 
the injected current. A proper analysis incorporates the Nernst and Seebeck effects considering temperature gradients 
in all three directions in space. We have treated these effects case-by-case to identify their symmetry and their possible 
contribution to RYJ ■ 
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The anomalous Nernst effect (ANE) is analogous to the anomalous Hall effect driven by a temperature gradient S5,S7 . 
The symmetry of the electric held driven by the ANE that gives rise to a second harmonic signal is as following: 


E AjVE = -aVTxM, (S12) 

where a is the material dependent ANE coefficient. By considering an arbitrary temperature gradient VT = 
{VT x ,VT v ,VT z ) and magnetization direction M = (sin 9 cos <p, sin 9 sin p, sin 9) we find the following allowed symme¬ 
tries for the second harmonic longitudinal {R%J) and transverse (i? 2 ^ VT ) resistances: 

R‘ 2 ^ E ck. X?T y cos 9 — VT Z sin 9 sin p , (S13) 

R 2 u ANE °c VT z sin0cosp- VT x cos6». (S14) 

The anisotropic magneto-thermopower (AMTEP) is the magnetization-dependent Seebeck effect and is the thermal 
analogous of the magnetoresistance S14 . Therefore, its manifestation follows Eq. S3 and, in terms of temperature 
gradients, reads: 


R 2 u ITEP c* VT x (sin 2 6* cos 2 p + sin 2 9 sin 2 ip ), (S15) 

R 2 oj AMTEP ^ VT !/ (sin 2 9 cos 2 ip + sin 2 9 sin 2 p). (S16) 

The transverse manifestation of the AMTEP, the so-called planar Nernst effect (PNE) S15,S16 , is analogous to the 
planar Hall effect and reads: 


R P ^ E ocVT y sin 2 9 sin 2ip + VXl sin 29 cos 2 p , (S17) 


R 2 u PNE « VT x sin 2 9 sin 2p + VT Z sin 29 sin 2 ip . (S18) 

In addition to these signals, in FM/NM bilayers one has to take into account the spin Seebeck effect (SSE), 
whereby a thermal gradient can drive a spin current that can be detected as a voltage across the bilayer via the 
inverse spin Hall effect. In our measurement geometry we could, in principle, observe the longitudinal SSE S9 due 
to a perpendicular thermal gradient, which would give rise to a signal with identical symmetry as that of the ANE 
(Eq. S13 and Eq. S14) S8 , although its sign would depend on the sign of the spin Hall effect in the NM layer. 

By considering Eq. S13 through Eq. S18 we notice that only a signal originating from the ANE (possibly including 
also the longitudinal SSE) and due to VT Z possesses the same symmetry as that of the USMR {R% u f MR oc sin0siny>). 
In the following, we discuss how these two effects can be separated. In previous work, we have shown that the ANE 
driven by a perpendicular temperature gradient induces a nonzero second harmonic Hall resistance in Co layers S6 . We 
use an 8 nm thick Co sample as a reference to reveal the influence of the ANE on the second harmonic longitudinal 
resistance. Figure S3 shows the first and second harmonic signals of longitudinal and transverse measurements when 
M is rotated in the xy plane in a fixed field B = 0.4 T. Sinusoidal fits (red solid curves) confirm the expected MR 
(planar Hall effect) symmetry in the longitudinal (transverse) first harmonic signals and the ANE symmetry in the 
second harmonic signals. Since the MR and planar Hall effect have the same microscopic origin, their amplitude is 
proportional to the physical distance over which they are measured. The same holds for the longitudinal and transverse 
ANE signals, which implies that the ratio between the longitudinal and transverse MR and ANE measurements is 
equal to the length/width (l/w) ratio of the Hall bar. This was verified experimentally by measuring the ratio l/w 
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FIG. S3. Longitudinal and transverse measurements of the MR and ANE in a single Co layer. First and second 
harmonic signals corresponding to the longitudinal and transverse resistance of Co(8nm). The sinusoidal fits (red curves) to 
the first harmonic signals show that the magnetization follows the applied field direction and the angular dependence of the 
MR (Ruj) and planar Hall effect whereas the second harmonic signals have the angular dependence expected of the ANE 

according to Eqs. S13 and S14. Note that — 4 = l/w as measured by scanning electron microscopy for this device. 



FIG. S4. Separation of AD-SOT and ANE contributions to L- Plot of the cosine component of L due to the AD-SOT 
and ANE as a function of ( B ex t + Bbem + Bpma )^ 1 measured in an xy scan for Ta(6nm)|Co(2.5nm) and Pt(6nm)|Co(2.5nm). 


of our Hall bars using scanning electron microscopy for devices of different size and including other reference layers 
such as Ti|Co bilayers for which the ANE and MR coefficients differ from those of a single Co layer S6 . Therefore, by 
measuring the thermal signal in the transverse geometry one can accurately determine its sign and magnitude in the 
longitudinal measurements. 

In FM/NM bilayers we must additionally consider the fact that the transverse thermoelectric signal, R^ u j VT , is 
mixed with R 2 U SOT ■ These two effects, however, can be separated in a quantitative way as shown in Ref. S6 and 
explained briefly below. During an xy scan of the magnetization, the FL-SOT gives a contribution to R^ proportional 
to cos 3 tp + cos tp, whereas the AD-SOT and the ANE both give a contribution proportional to cos <p. The AD-SOT 
and ANE contributions can be further separated by considering that the AD-SOT induces dynamical oscillations 
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of the magnetization, the amplitude of which is proportional to the magnetic susceptibility of the FM layer. The 
resulting OT signal therefore depends on the susceptibility of the magnetization, which decreases with increasing 
external field, whereas the ANE contribution is constant provided that the magnetization is saturated along the field 
direction, as is the case in our measurements. The cos <p component of is i n fact a linear function of the inverse of 
the effective magnetic fields acting on the magnetization, ( B ext + Bdem + Bpma ) , with slope proportional to the 

AD-SOT and intercept proportional to the ANE S6 . Figure S4 shows a plot of this component for two of the samples 
used in this study, Ta(6nm)|Co(2.5nm) and Pt(6nm)|Co(2.5nm). The intercept of the linear fit gives a non-negligible 
ANE-induced signal in the case of Ta {R^ u j VT = 11 mil), whereas a negligible ANE is found for the Pt sample 
= 0.02 mil). The corresponding AIR of thermal origin is then given by R%J = ;g7?2w VT - The above study is 
repeated for all the samples used in this work to separate the thermoelectric contributions from the measured second 
harmonic MR and reveal the pure USMR signal, given by R 2 RMR = — R%J ■ 


SI 5. ABSENCE OF THE USMR IN SINGLE LAYER SAMPLES 


The measurements presented in Fig. S3 exclude the existence of a net USMR signal for a single Co layer. In order 
to further verify if the USMR is due to the Co/NM interface or to another nonlinear mechanism giving rise, e.g., 
to a self-induced MR induced by the SHE, such as that due to edge spin accumulation predicted by Dyakonov S17 , 
we have grown two Ta(6nm) and Pt(6nm) reference layers and measured R u and for these samples. Figure S5 
shows the zy scans for Ta (upper panels) and Pt (lower panels). The measurements are performed using a current 
density j = 0.2 • 10' A/cm 2 for Ta and j = 10 7 A/cm 2 for Pt. The lower current density used for the Ta layer is 
due to the input voltage limit of our acquisition system and the high resistivity of Ta. Within the sensitivity of our 
measurements, we do not observe any systematic magnetic field dependence of either R u or i? 2 ui> which excludes any 
influence from the NM layer alone in the USMR. 
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FIG. S5. Measurements of R u or in single Ta and Pt layers. First (left) and second (right) harmonic longitudinal 
resistance signals in Ta(6nm) (top) and Pt(6nm) (bottom) reference layers. Due to the high resistivity of Ta, the Ta data are 
averaged over five zy scans performed with B ex t = 1.7 T and j = 0.2 • 10' A/cm 2 . For Pt, a single zy scan is shown with 
j = 10 7 A/cm 2 and B ex t = 1.9 T. 















SI 6. USMR IN Ta|Cu|Co LAYERS 


In order to elucidate the possible role of the FM-induced magnetization in the NM, we have inserted a Cu spacer 
with thickness 2 and 4 nm between Ta and Co. The choice of Cu is motivated by its large spin diffusion length as 
well as by its low and short-range induced magnetization. Ta(6nm)|Co(2.5nm) was chosen as a test system because 
of the larger absolute amplitude of the USMR relative to Pt. The aim of this study is twofold. First, to see if the 
USMR is still present when the Ta and Co layers are physically separated. Second, to examine if there is any influence 
of the induced magnetization on the reported USMR signal, although, judging by the change of sign of the USMR 
observed in the Ta|Co and Pt|Co layers, this appears unlikely. If the induced magnetization plays a role, we do 
not expect to measure a USMR signal after the insertion of the Cu spacer layers. However, if the USMR is due to 
the current-induced spin accumulation at the FM/NM interface, we expect to measure R 2 RMR 7^ 0 but reduced in 
amplitude with respect to Ta|Co, mainly due to current shunting in the Cu layer. 



cp B (degrees) 

FIG. S6. Influence of a Cu spacer layer on the USMR signal of Ta|Co. xy scan of R 2 RMR measured in Ta(l,6nm)|Cu(2- 
6nm)|Co(2.5nm) trilayers with current density j = 10 7 A/cm 2 . 

Figure S6 shows the USMR of Ta(l,6nm)|Cu(2-6nm)|Co(2.5nm) layers. We observe a clear signal in Ta(6nm)| 
Cu(2,4nm)|Co(2.5nm) that has the same sign as that measured in Ta(6nm)|Co(2.5nm). As expected, the USMR 
decreases significantly with increasing Cu layer thickness, from 1.9 ml! for Ta(6mn)|Cu(2nm)|Co(2.5nm) to 0.28 mil 
for Ta(6nm)|Cu(4nm)|Co(2.5nm), and so does the resistivity of the devices, from 54.6 /iHem to 34.1 pflcm, showing 
that the shunting effect of Cu quickly dominates the conduction. The sample with Ta(lnm) has a resistivity of 
16.9 ^flcm and serves as a reference to show that the USMR vanishes if the conduction is dominated by the Co/Cu 
bilayer, due to the negligible SHE of Cu. We note that the 1 nm of Ta was deposited as an adhesion layer for Cu on 
the SiO substrate and that its conductivity is expected to be negligibly small due to oxidation and the much larger 
resistivity of Ta relative to Cu. Overall, these data confirm that the USMR is dominated by the spin accumulation at 
the FM/NM interface induced by the SHE in the heavy metal layer and that the induced magnetization in the NM 
does not play a significant role. 


SI 7. SPIN ACCUMULATION AND POTENTIAL SHIFT INDUCED BY THE SHE IN A FM/NM 

BILAYER 

There are essentially two effects that are important to interpret the USMR and model its dependence on the 
thickness of the NM layer. One is the buildup of the SHE-induced spin accumulation at the FM/NM interface, which 
occurs over a length-scale comparable to the spin diffusion length of the NM, and the other is a dilution effect due 
to reduced number of electrons that scatter at the interface relative to the total number of electrons participating to 
the conduction. We consider these two effects separately. 
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Our first point is the assumption that the USMR is proportional to the nonequilibrium spin accumulation at 
the FM/NM interface induced by the in-plane charge current. As explained in the text, this is motivated by the 
modulation of the conductivity mismatch that exists at the interface between a FM and a NM due to the SHE- 
induced change of majority and minority spin populations. In order to calculate the spin accumulation, we adopt the 
drift-diffusion approach 818,819 . The spin accumulation is defined as / j, s = — where /x and //A represent the spin- 

dependent electrochemical potentials for majority (f) and minority (j.) electrons. We recall that the nonequilibrium 
magnetization in the NM is 5m = — N(£p)yLBP s i where A !'{sf) is the density of states at the Fermi level and the 
minus sign stems from the opposite orientation of the magnetic spin moment and spin angular moment. 

We consider a bilayer consisting of a NM with electrical conductivity apr, spin diffusion length Ajv, and thickness 
t n , and a FM defined by the analogous quantities a F , X f , and t F . The conductivity of the FM is assumed to be 
the sum of the independent conductivities for the majority and minority electrons: <jf = cA + crk Accordingly, the 
current flowing in the FM has a net spin polarization P = (erf — cr^)/(cA + erf). We define x as the current direction 
and ^ as the direction normal to the interface; the interface plane is situated at z = 0. The source term for /x s is the 
spin current generated inside the NM by the SHE and propagating along 2 : 

j° = 3sh (J x z) = ~j ° SH y = - 6 S H<TNE x y , (S19) 

where y represents the spin polarization direction, E x the electric field driving the charge current j through the bilayer, 
and Osh the spin Hall angle of the NM. The spin Hall angle of the FM is assumed to be zero. This approximation 
neglects the spin and charge accumulation induced by the anomalous Hall effect at the boundaries of the FM, which 
can be included in more detailed calculations but are not essential to the arguments developed here. Furthermore, 
we assume that the magnetization of the FM is saturated parallel to ±y. In this case, there are no sink terms for the 
spin current inside the FM other than spin-flip relaxation; in particular, the spin current associated to the SOT via 
the real and imaginary part of the spin-mixing conductance (G-jq) vanislres S2 °: 

jf OT - Re{G t j.}r J 4 D + Im{G n }r F L ~ Re{G t4 .}M x (y x M) + Im{G n }(M x y) = 0. (S20) 


With this assumption the spin diffusion equation governing the spin accumulation reduces to a one-dimensional 
problem: 


^ Msf ,n ( z ) 


Msf.jv ( Z ) 
^ 2 F,N 


(S21) 


where Hs f ,n is the y component of the spin accumulation vector in either the FM 
j s = js y is then given by 

Jsn 0 ) = - /4r) - Jsh in the NM 

js F (z) = -\ d z Af - in the FM 

The general solutions of Eq. S21 read 

y] v (z) = A N + B n z + C N exp(z/X N ) + D n exp(-z/X N ) 
An( z ) ~ A’n + B' n z + C' N exp( 2 / Xm) + D' N exp(— zjX n) 
nlr(z) = A f + B F z + Cp exp(z/Xp) + Dp exp(— z/Xp) 
Hf(z) = A'p + B'pZ + C'p exp(z/X F ) + D' F ex^,(-z/X F ) 


or the NM layer. 

The spin current 

layer (z < 0), 

(S22) 

layer (z > 0). 

(S23) 

Jz < 0 

(S24) 

jz > 0 

(S25) 


The coefhcients appearing in Eqs. S24, S25 are determined by imposing charge conservation V • j = 0 in the two layers: 

= ° z<0 ’ ( S26 ) 

d^(a^/j^p + crf Af) =0 z > 0, (S27) 


and by the following boundary conditions: 
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1. The net charge flow perpendicular to the FM/NM interface is zero: 

3f + 3f = 3n + 3n = ( S28 ) 

2. The spin-dependent current components are continuous across the FM/NM interface: 

J^(0)=jV(0). (S29) 

3. The discontinuity of the electrochemical potentials of majority and minority electrons at the FM/NM interface 
is proportional to the spin-dependent boundary resistances b and Hq respectively: 

Mf( 0) - Aijv(°) = er T j T (0) and p l F (0) - p l N (0) = er l j l ( 0). (S30) 


4. The spin current vanishes at the bottom of the NM layer and at the top of the FM layer: 

js N {~tN ) = js F (t F ) = 0. 


(S31) 


In Eqs. S28-S30 above we have introduced the majority and minority current components, j p = / d z p,' p and 

jjv*' = =F \j%H> where the ’’majority” direction in the NM is defined with respect to that of the FM. The 

boundary condition S29 implies the continuity of the spin current, j SF (0) = j SN (0), which is justified only in the 
absence of spin-flip scattering at the interface 319 . To simplify our model, we assume this to be the case even if spin- 
flip scattering at the Co|Pt interface can reduce the spin current by as much as a factor two S21 . In Eq. S30 we retain 
the effects of interface scattering due to the different band structure of the NM and FM and diffuse scattering by the 
roughness or chemical disorder of the interface, which are relevant in the theory of CIP-GMR S22 . The spin-dependent 
boundary resistances are defined in analogy with the spin-dependent resistivity in bulk FM: 


r t(l) = 2 r b( l - (+b)> 


(S32) 


where rf, = (r-j- + rj .)/4 is the interface resistance (in units of flm 2 ) and 7 = (77. — G")/( r t + r t) the interfacial spin 
asymmetry coefficient 319 . 

Finally, we obtain 


AGn ( 0 ) = tanh 


t N 1 + ^( 1 - J32 )tanh^ 

2A W ! + CQth ^ (! - P 2 ) tanh £ ’ 


(S33) 


where pP s ^, = 2 e pn^n Osh j is the bare spin accumulation due to the SHE that would be obtained for a single, infinitely 
thick NM layer and we have substituted the conductivity by the corresponding resistivity parameters pn = 1/<r n and 
pF = Note that in the limit pF —> 00 we recover the spin accumulation calculated for an insulating FM/NM 

interface 320,323 . Differently from the latter case, however, the SHE induces a shift of the electrochemical potential of 
the NM layer, p N = (p? N + p^)/ 2 , relative to that of the FM, p F = (p} p + p p )/ 2 . Taking p F ( 00) = 0 as the reference 
level, we have 


Ap N = p N ~ Pf{ 00 ) = ~(P + lf)p° tanh 


tN 1 


1 


2A tv 1 T 


1 + 


PN * N (1 — P 2 ) tanh 1/’. coth 4^ 
PF X F __ A y _ A jv 

1 — r 


(S34) 


where f = p //( F (1 — P 2 ) tanh 

By measuring the resistivity of Co, Pt, and Ta reference layers and using published values for P, Atv, A p, 
and 7 , it is possible to estimate the relative weight of the different terms appearing in Eqs. S33 and S34. For the 
resistivity, we obtain pc 0 = 36.1 /rflcm, pp t = 33.4 pttcm, and pp a = 228 pflcm. We take P = 0.31 and Ac*o — 30 nm 
from CPP-GMR measurements of Co|Cu films with resistivity similar to ours 324 , rj, = 0.74 film 2 and 7 = 0.53 as 
recently measured for Co|Pt (Ref. S21). Lacking similar measurements for Co|Ta, we assume the same interface 
resistance values as for Co|Pt. Further, we take A Pt « \ Ta « 1.5 nm (Refs. S25 and S26). With these parameters 
and tp = 2.5 nm, we calculate the products (1 — P 2 ) tanh = 0.0034 (0.023) for Co|Pt (Co|Ta) and r = 0.005. 
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FIG. S7. Spin accumulation and modulation of the spin-dependent electrochemical potential by the SHE. Profile 
of the electrochemical potential of majority (pJ, blue lines) and minority (/A, red lines) electrons in proximity of the FM/NM 
interface calculated using Eqs. S24-S25 for a, positive and b, negative current. The electrochemical potential of the NM shifts 
relative to that of the FM as indicated by the green arrow. The direction of the magnetization is M || y in both panels. 


Inserting these products into Eqs. S33 and S34, we observe that, for fjv <1 Aw, the spin accumulation and the 
electrochemical potential shift of Co|Pt (Co|Ta) are approximated to within 1% (3%) by 

/W(0) ~ ^°~ tanh 2l^> ( S35 ) 

and 

A p N ~ -Pp° a „ tanh ( S36 ) 

More generally, we notice that, as long as tpr > Aw, tp < A p, Pn^n < PfAf, and rb < PfXf, both p SN ( 0) and A pjy 
are mainly determined by the properties of the NM. Figure S7 illustrates the behavior of p\ F , F , and A p, N as 
a function of 2 calculated using Eqs. S24-S25 and the following parameters: j = 10 7 A/cm 2 , d SH = 0.1, P = 0.5, 
Pn = Pf = 30 /iflcm, Aw = 1.5 nm, A p = 10 nm, Tb = 7 = 0. Figure 5 of the main text reports a similar calculation 

where rb = 0.74 fflm 2 and 7 = 0.53 as appropriate for Co|Pt (Ref. S21), which introduces a discontinuity between 

/zjy' and Pp^ at z = 0. 


SI 8. CURRENT SHUNTING THROUGH THE NM AND FM LAYERS 

As the USMR is a pure interface effect, we expect a decrease of /±r usmr when the thickness of the NM and FM 
layers becomes larger than the respective spin diffusion lengths and the current is shunted away from the interface. 
This effect is similar to the ’’dilution” of the CIP-GMR observed in FM/NM/FM multilayers 827 . However, differently 
from CIP-GMR, the relevant lengthscale beyond which the dilution effect becomes significant here is the spin diffusion 
length rather than the electron mean free path. 

We use a simple parallel resistor model to describe the current flow in a FM/NM bilayer. In this model, the 
’’inactive” regions of the NM and FM layers are described by the resistances Rn and R F whereas the ’’active” interface 
region is described by a linear, current-independent resistance Ri in series with a nonlinear, current-dependent USMR 
resistance r s . The sign of r s is determined by the cross product j x M. Retaining only first order terms proportional 
















12 


to r s , the equivalent resistance of Rn, Rp, and Rj ± r s connected in parallel is 

Ri 


R ± = 


1 + ^^ E±rS ’ 


which gives 


A R! 


USMR 


= R + - Rr = 


2 r s 


( 1 + ^^ E ) 


2 ' 


(S37) 


(S38) 


To analyze the USMR dependence on the thickness of the NM layer alone, it is practical to include the constant 
resistance terms Rp and Ri in a single parameter Rpj = RpRi/{Rp + Rj). We thus obtain 


A R 


USMR 


2 r. 


(l+Tff tN 


2 ’ 


(S39) 


and 


where 


ARUSMR 2 r. 

R 


Rfi 1 + AAA T tN ' 


R = 


R F + R RnRfi 


Rfi 


2 R n + R fi l + R£J-ft N ’ 

Pn t 

and w/l is the ratio between the width and the length of the current line. 


(S40) 


(S41) 


SI 9. DEPENDENCE OF THE USMR ON THE THICKNESS OF THE NM LAYER 


As argumented in the main text and above, we expect the USMR to be proportional to A/i times a dilution factor 
accounting for the shunting of the charge current by the NM layer. Combining Eqs. S35, S39, and S40 we obtain the 
following phenomenological expressions: 


A R 


USMR. 


= A- 


tanh JA- 

_ 

f l + teLV t N 

\ Pn i 


2 ’ 


ARU smr 
R 


A 


tanh 7pr- 

ZAn 


RFll + ARTT t *' 


(S42) 


where A is a fitting parameter proportional to Ppn^nOshJ that accounts also for quantitative differences entering 
into the calculation of the resistance that escape our model. The values of the free parameters of the fits presented in 
Fig. 4 of the main text are A = 7 mil (-14 mil) for Ta|Co (Pt|Co), Rpi/ppr = 0.08 nm -1 (0.44 mn" 1 2 3 4 5 6 ), X Fa = 1.4 nm 
and Xpt = 1.1 nm. 
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